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We present an improved measurement of CP-violating asymmetries in B° — > tt + it~ decays 
based on a 78 fb _1 data sample collected at the T(AS) resonance with the Belle detector 
at the KEKB asymmetric-energy e + e~ collider. We reconstruct one neutral B meson as 
a B° — > tt + tt~ CP eigenstate and identify the flavor of the accompanying B meson from 
its decay products. From the asymmetry in the distribution of the time intervals between 
the two B meson decay points, we obtain the CP-violating asymmetry amplitudes Ann = 
+0.77 ±0.27(stat) ±0.08(syst) and Snn = -1.23 ±0.41(stat) ±S;JJ?(syst), where the statistical 
uncertainties are determined from Monte Carlo pseudo-experiments. We rule out the CP- 
conserving case, Ann = Snn = 0, at the 99.93% confidence level. We discuss how these results 
constrain the value of the CKM angle cj>2- 



1 INTRODUCTION 

Kobayashi and Maskawa (KM) proposed, in 1973, a model where CP violation is incorporated 
as an irreducible complex phase in the weak-interaction quark mixing matrix^. The KM model 
predicts CP-violating asymmetries in the time-dependent rates for B° and B° decays to a 
common CP eigenstate, fcp, as 

A C p{t) = T(W - fcp . t) - T(B0 - — t) = Af cos(Am^) + S f sm(A^t), (1) 

where T(B° f C p;t)(T(B° fcp\t)) is the partial decay rate of pure B°(B°) state at t = 0, 
and Am,} is the mass difference between the two B° mass eigenstates ™. The CP- violating 
parameters Af and Sf defined in Eq. Q are expressed by Af = (\Xf\ 2 — 1)/(|A/| 2 + 1) and 
Sf = 2Im\f / ( | Ay j 2 + 1), where Aj is a complex parameter that depends on B°B° mixing, 



on CP eigenvalue of fcp and on the amplitudes for B° and B° decay to fcp- In the mode 
B° — > J/ipKs ^, the Standard Model predicts Sj/$k s = sin20i - and Aj/^j^ = 0. Recent 
measurements of the CP- violating parameter sin20i by the Belled and BaBar collaborations 
established CP violation in the b — > ccs that is consistent with KM expectations. In the mode 
B° — > 7r + 7r _ , we would have S n7T = sin 202 and A^n = 0, or equivalently |A/| = 1, if the b — > u 
tree amplitude were dominant. The situation is complicated by the possibility of significant 
contributions from gluonic b — > d penguin amplitudes that have a different weak phase and 
additional strong phases 1 -^. As a result, may not be equal to sin 2cp2 and direct CP violation, 
■A.-K1T 7^ 0, may occur. 



2 EXPERIMENTAL APPARATUS 

In this paper, we report an updated measurement H that is based on a 78 fb 1 data sample, 
containing 85xl0 6 BB pairs, w hich were collected with the Belle detector ^ at the KEKB 
asymmetric-energy e + e~ collider^^ at T(4S) resonance. KEKB archives the highest peak lu- 
minosity of 8.26 x 10 33 cm -2 s -1 as of March, 2003. The Belle detector is a large-solid-angle 
general purpose spectrometer that consists of a silicon vertex detector (SVD), a central drift 
chamber (CDC), an array of aerogel threshold Cerenkov counters (ACC), time-of-flight scintil- 
lation counters, and an electromagnetic calorimeter comprised of CsI(Tl) crystals located inside 
a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux return located 
outside of the coil is instrumented to detect K® mesons and muons. 



3 EXPERIMENTAL METHOD 

Experimental method in this report is very similar to the well-established sin 2^»i measurement^ 
in B° — » ccKs- In the decay chain T(45) — > B°B° — > fcpft&g, where one of B mesons decays 
at time tcp to fcp and the other decays at time tt a g to a final state /tag that distinguishes 
between B° and B°, the decay rate has a time dependence given by^^ 

e -|At|/T B o 

V^(At) = — [1 + q ■ {S n7T sin(Am d At) + A™ cos(Am d At)}] , (2) 

4t b o 

where t b o is the B° lifetime, At = tcp — ttagi an d the 6-flavor charge q = +1 (—1) when the 

tagging B meson is a B (B°). As described below in Section |IJ we reconstruct one of B mesons 

with two tracks identified as pion. We determine the flavor of the accompanying B from the 

information of its decay products. At is determined from Az, the displacement in z between 

the decay vertices of two B mesons: At ~ [zqp — zt £Lg )/(3'yc = Az/fi'yc, where is a Lorentz 

boost factor of 0.425. In flavor t aggi ng and At measurement, we apply the same method used 

for the Belle sin20i measurement^^. 

The CP asymmetry parameters, A nn and S n7T , are obtained from an unbinned maximum 

likelihood fit to the observed At distribution. For this purpose, we use probability density 

functions, Pj (PDFs) that are based on theoretical distributions given by Eq.(|2"|). PDFs are 

diluted with the background and smeared by the detector response. We use the same detector 
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response as those used for the sin20i measurement 1 1 • In the fit, S n7r and A^n are free 
parameters determined by maximizing the likelihood function C = YliPi, where the product is 
over all B° — > ir + ir~ candidates. 



4 B° -> vr+vr RECONSTRUCTION 

The B° — > 7r + 7r _ event selection is described in detail elsewhere^! We use oppositely charged 
track pairs that are positively identified as pions according to the combined information from the 



ACC and the CDC dE/dx measurement. The efficiency of the pion identification is 91% and the 
contamination from miss-identified kaon is 10% for the track momentum from 1.5 to 4.5 GeV/c. 
Candidate B mesons are reconstructed using the energy difference AE = E c B ms — E^ CSLm CTas 

and the beam-energy constrained mass M hc = \J {E£ cam cms ) 2 - (p§" s ) 2 , where E£ eam cms is the 
cms beam energy, and E™ 8 and p§° s are the cms energy and momentum of the B candidate. 
The major sources of the background in B° — > 7r + 7r~ are B° — > K + ir~ decays where kaons are 
misidentified as pions, and the background from the e + e~ — > qq continuum (q = u, d, s, c). 

In order to suppress the continuum background, we form signal and background likelihood 
functions, £s and £bg, from t wo v ariables. One is a Fisher discriminant determined from six 
modified Fox- Wolfram moments'-^, the other is the cms B flight direction wi th r espect to the z 
axis. We determine £s from a GEANT-based Monte Carlo (MC) simulation^!, and £bg from 
AE-Mb c sideband data dominated by the continuum background. We reduce the continuum 
background by imposing requirements on the likelihood ratio LR = £s/{£s+£bg) f° r candidate 
events. We optimize LR requirement for each flavor tagging category to maximize the expected 
sensitivity. 

5 SIGNAL YIELD 

Figures ^a) and (b) show the AE distributions for the B° — > ir + ir~ candidates that are in the 
Mb c signal region with LR > 0.825 and with LR < 0.825, respectively, after flavor tagging 
and vertex reconstruction. In the Mt, c and AE signal region, we find 275 candidates for LR 
> 0.825 and 485 candidates for LR < 0.825. The B° -> vr+Tr" signal yield for LR > 0.825 
is extracted by fitting the AE distribution with a Gaussian signal function plus contributions 
from misidentified B° — > K + ir~ events, three-body .B-decays, and continuum background. The 
fit yields 106^5 ir + ir~ events, 4llg° K + tt~ events and 128]t 6 continuum events in the signal 
region. The errors do not include systematic uncertainties unless otherwise stated. Here the 
error on the yield of continuum events in the signal region is obtained by scaling the error of 
the yield from the fit that encompasses the entire AE range. For LR < 0.825, we fix the level 
of 7r + 7r - signal by scaling the LR > 0.825 number by a MC-determined factor and that of the 
continuum background from the sideband. The ratio of the K + tt~ background to the ir + ir~ 
signal is fixed to the value measured with the LR > 0.825 sample. We obtain 57 ± 8 ir + ir~ 
events, 22^ K + tt~ events and 406 ± 17 continuum events in the signal region for LR < 0.825. 
The contribution from three-body i?-decays is negligibly small in the signal region. 

6 FIT RESULTS 

The unbinned maximum likelihood fit to the 760 B° — > ir + ir~ - candidates (391 B°- and 369 
B°-tags), containing 163l 2 3 7r+7r_ signal events, yields 

= +0.77 ±0.27(stat) ±0.08(syst), (3) 
S n7T = -1.23 ±0.41(stat) t^(syst). (4) 

Here we quote the rms values of the MC and S n7T distributions as the statistical errors of 
our measurement, because the log-likelihood ratio curves from our data deviates from parabolic 
behavior and the statistical error estimation using log-likelihood ratio, — 21n(£/£ max ), is not 
appropriate. 

7 CROSS CHECKS 



We perform a number of cross checks. We measure the B meson lifetime using the same vertex 
reconstruction method. In addition, we check for biases in the analysis using samples of non-CP 



eigenstates, B° — > K + ir~ decays, and sideband data. 

We perform a B° lifetime measurement with the B° — > tt + it~ candidate events that uses the 
same background fractions, vertex reconstruction methods, and resolution functions that are 
used for the CP fit. The result, r B o = 1A2±^ ps, is consistent with the world-average valueP. 
We also perform measurement of B° lifetime and B°-B° mixing with the 1371 B° — > K + ir~ 
candidate events (610 signal events), that are selected using positively identified the charged 
kaons. we use the same vertex reconstruction method and wrong-tag fractions as B° — > 7r + 7r~ 
and determine t b o = 1.46 ±0.08 ps and Am^ = 0.55t_Q^ ps -1 , which are in agreement with the 
world average values^! 

We also use samples of non-CP eigenstate B° — > D~ir + , D*~ir + and D*~ p + decays, selected 
with the same event-shape criteria, to check for biases in the analysis. The combined fit to 
this control sample of 15321 events yields A = -0.015 ± 0.022 and S = 0.045 ± 0.033. As 
expected, neither mixing-induced nor direct CP-violating asymmetry is observed. A fit to the 
B° — » K + ir~ candidates candidates yields Ak-k = —0.03 ±0.11, in agreement with the counting 
analysis mentioned above^J and Sk-k — 0.08±0.16, which is consistent with zero. A comparison 
of the event yields and At distributions for B°- and P°-tagged events in the sideband region 
reveals no significant asymmetry. 



8 SIGNIFICANCE 

We use the Feldman-Cousins frequentist approach^! to determine the statistical significance of 
our measurement. In order to form confidence intervals, we use the A n7T and distributions 
of the results of fits to MC pseudo-experiments for various input values of A^n and S nn . The 
distributions incorporate possible biases at the boundary of the physical region as well as a 
correlation between A n7T and S^] these effects are taken into account by this method. The dis- 
tributions are also smeared with Gaussian functions that account for systematic errors. Figure El 
shows the resulting two-dimensional confidence regions in the A^-n vs. S n7T plane. The case that 
CP symmetry is conserved, A^ = S W7T = 0, is ruled out at the 99.93% confidence level (C.L.), 
equivalent to 3.4a significance for Gaussian errors. The minimum confidence level for A n7T = 
0, the case of no direct CP violation, occurs at (<S W7r , A n7T ) = (—1.0,0.0) and is 97.3%, which 
corresponds to 2.2a significance. 



9 DISCUSSION 

Using the standard definitions of weak phases <j)\, 02, and 03, the decay amplitudes for B° and 
B to ir + Tr~ are 

A{B°^n + 7r-) = -(\T\e iST e i(t>3 +\P\e i5p ), (5) 
^(P ^vr + vr-) = -(\T\e i&T e- i(t>3 + \P\e iSp ), (6) 

where T and P are the amplitudes for the tree and penguin graphs and St and 5p are their 
strong phases. Here we adopt the notation of and use the convention in which the top- 

quark contributions are integrated out in the short-distance effective Hamiltonian. In addition, 
the unitarity relation V* b V u d + V*b^cd = —Vtb^td is applied. Using the above expressions and 
02 = vr - 0i - 3 , we determine = e 2 ^ 2 [l + \P/T\e i{ - s+ ^]/\l + \P/T\e^ s ^}. Explicit 



expressions for <S T7r and A nn are 

S n7T = [sin20 2 + 2|P/T|sin(0i - 2 )cos5 - |P/T| 2 sin20i]/"K, (7) 

Attt = -[2[P/r|sin(0 2 ±0i)sin5]/^, (8) 

K = l-2|P/T|cos<5cos(0 2 ±0i)±|P/T| 2 , (9) 



where 5 = 5 P - 5 T . We take -180° < 5 < 180°. When A^ is positive and 0° < fa + fa> < .180° 5 
is negative. Recent theoretical estimates prefer \P/T\ ~ 0.3 with large uncertainties^! '^ 
Figures Eta)-(e) show the regions for </>2 and 8 corresponding to the 68.3% C.L., 95.5% C.L. and 
99.73% C.L. region of A^-n and (shown in Fig. EJ) for representative values of \P/T\ and 
Note that a value of (S n7T ,A n7T ) inside the 68.3% C.L. contour requires a value of \P/T\ 
greater than ~0.3. The allowed region is not very sensitive to variations of fa within the errors 
of the measurements, as can be seen by comparing Figs.^a), (c) and (e). The range of 4>2 that 
corresponds to the 95.5% C.L. region of A^ and S n7T in Fig. 01 is 

78° < cp 2 < 152°, (10) 

for fa = 23.5° and 0.15 < \P/T\ < 0.45. Th e result is in agreement with constraints on the 
unitarity triangle from other measurements^^. 



10 CONCLUSION 



In summary, we have performed an improved measurement of CP violation parameters in B — > 
7r + 7r _ decays. An unbinned maximum likelihood fit to 760 B° — > 7r + 7r~ candidates, which 
contain 163^23 (stat) ir + ir~ signal events, yields A^k = +0.77 ± 0.27(stat) ± 0.08(syst), and 
S niT = —1.23 ±0.41(stat) ^QQy(syst), where the statistical uncertainties are determined from MC 
pseudo-experiments. This result is consistent with our previous measurement^^ and supersedes 
it. We obtain confidence intervals for CP-violating asymmetry parameters A nrr and based on 
the Feldman-Cousins approach where we use MC pseudo-experiments to determine acceptance 
regions. We rule out the CP-conserving case, A^n = = 0, at the 99.93% confidence level. 

The result for S nn indicates that mixing-induced CP violation is large, and the large A^ 
term is an indication of direct CP violation in B meson decay. Constraints within the Standard 
Model on the CKM angle <p2 and the hadronic phase difference between the tree (T) and penguin 
(P) amplitudes are obtained for \P/T\ values that are favored theoretically. We find an allowed 
region of <p2 that is consistent with constraints on the unitarity triangle from other measurements. 
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Figure 1: AE distributions in the Mt, c signal region for (a) B° — > 7r + 7r~ candidates with LR > 0.825 and (b) 
B° — > 7r + 7r~ candidates with LJ? < 0.825. The sum of the signal and background functions is shown as a solid 
curve. The solid curve with hatched area represents the tt + tv~ component, the dashed curve represents the K + n~ 
component, the dotted curve represents the continuum background, and the dot-dashed curve represents the 

charmless three-body B decay background component. 




Figure 2: The raw, unweighted At distributions for the 275 B — > tv 7v~ candidates with LR > 0.825 in the 
signal region: (a) 148 candidates with q = +1, i.e. the tag side is identified as B°; (b) 127 candidates with q = — 1; 
(c) B° — > 7v + tt~ yields after background subtraction. The errors are statistical only and do not include the error 
on the background subtraction; (d) the CP asymmetry for B° — > ir + n~ after background subtraction. In Figs, 
(a) through (c), the curves show the results of the unbinned maximum likelihood fit to the At distributions of the 
760 B° — > tv + tv~ candidates. In Fig. (d), the solid curve shows the resultant CP asymmetry, while the dashed 
(dotted) curve is the contribution from the cosine (sine) term. 




Figure 4: The regions for fa and 5 corresponding to the 68.3%, 95.5%, and 99.73% C.L. regions of Ann and S n n 
in Fig. for (a) fa = 25.9°, |P/T|=0.3, (b) fa = 23.5°, |P/T|=0.15, (c) fa = 23.5°, |P/T|=0.3, (d) fa = 23.5°, 
|P/T|=0.45, and (e) fa = 21.3°, \P/T\=0.3. The horizontal dashed lines correspond to fa = 180° - fa. 



